Homoclinic Orbits around Spinning Black Holes I: 
Exact Solution for the Kerr Separatrix 
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Under the dissipative effects of gravitational radiation, black hole binaries will transition from 
an inspiral to a plunge. The separatrix between bound and plunging orbits features prominently 
in the transition. For equatorial Kerr orbits, we show that the separatrix is a homoclinic orbit 
in one-to-one correspondence with an energetically-bound, unstable circular orbit. After providing 
a definition of homoclinic orbits, we exploit their correspondence with circular orbits and derive 
exact solutions for them. This paper focuses on homoclinic behavior in physical space, while in a 
companion paper we paint the complementary phase space portrait. The exact results for the Kerr 
separatrix could be useful for analytic or numerical studies of the transition from inspiral to plunge. 



I. INTRODUCTION 



A. Background and Motivation 

A direct observational detection of gravitational waves 
- perhaps the most fundamental prediction of a theory 
of curved spacetime - looms close at hand. Stellar mass 
compact objects spiraling into supermassive black holes 
have received particular attention as sources ofgravita- 
tional radiation for the planned LISA mission [lj. A di- 
rect detection of these extreme mass ratio inspirals (EM- 
RIs), as well as extraction of astrophysics @, y, HI, IE Is 13 > 
requires a thorough knowledge of the underlying dy- 
namics; it is the motion of the two bodies that shapes 
the gravitational waveform. A well-established approach 
models the EMRI as an adiabatic progression through a 
series of Kerr geodesies 0, 0, IE IE S EH ■ A transparent 
depiction of geodesic motion around spinning black holes 
is therefore essential, yet seemingly complicated [IE [H| 
and benefits from crucial signposts in the orbital dynam- 
ics. 

We decipher such a crucial signpost here. In partic- 
ular, we discuss an important family of separatrices in 
Kerr dynamics: the homoclinic orbits. 1 Around black 
holes, the homoclinic orbits are those that asymptotically 
approach the same unstable circular orbit in both the in- 
finite future and the infinite past, 2 as shown on the right 
of Fig. [TJ Under the identifier "separatrix" , homoclinic 
orbits have already garnered attention in the black hole 
literature [13 , Il3j - the homoclinic orbit is the separatrix 
between orbits that plunge to the horizon and those that 
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1 The terms "homoclinic orbit" and "separatrix" are, in this con- 
text, entirely interchangeable, although the former finds more 
use in the dynamical systems literature and the latter in the 
black hole and gravitational wave literature. 

2 Orbits that approach two different orbits in the infinite future 
and past, in contrast, are called heteroclinic orbits. 





FIG. 1: Left: A zoom-whirl orbit. Right: A homoclinic orbit 
approaching an unstable circular orbit. 



do not. The scenario of quasi-circular inspiral through 
a last stable circular orbit is a special example of the 
transition through a zero eccentricity homoclinic orbit. 
Orbits that merge before they have a chance to circular- 
ize will transit through an eccentric homoclinic orbit of 
the underlying conservative dynamics. Any analysis of 
the transition from inspiral to plunge will thus run into 
this special family. 

Homoclinic orbits are also a significant signpost for 
zoom-whirl behavior; an extreme form of perihelion pre- 
cession wherein trajectories zoom out into quasi-elliptical 
leaves en route to apastron and then execute multiple 
quasi-circular whirls near periastron before zooming out 
again, as shown in the left panel of Fig. [1] Though zoom- 
whirl behavior is sometimes thought to be associated only 
with highly eccentric orbits near the separatrix, we de- 
veloped a topological criterion for whirliness in and 
showed that in the strong-field regime orbits of any eccen- 
tricity can exhibit zoom- whirliness. Indeed, zoom-whirl 
behavior is neither exotic nor rare in the strong field [HI] • 
Still, homoclinic orbits are relevant as an infinite whirl 
limit in the distribution of geodesies, a connection we 
forge in this paper. 

Homoclinic orbits are therefore significant in shaping 
the geography of black hole orbits. In this first paper in 
a two part series, we devote some labor to resolving this 
landmark in physical space for equatorial orbits. (We 
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leave for a future work the non-equatorial case.) The 
pinnacle is an exact solution for equatorial homoclinic 
trajectories. A rarity among relativistic orbits, the exact 
solution can make semi-analytic treatment of the eccen- 
tric transition to plunge more wieldy. In paper II [Til ], 
we describe the flipside of the coin and detail the phase 
space portrait of the homoclinic orbits. We hope the re- 
sults will provide cohesion to the dynamical conversation. 

We begin this paper by finding exact expressions for 
the orbital parameters of the separatrices and use them 
to derive Eqs. (|26| . exact expressions for the trajectories 
themselves. As this paper is concerned with the physical 
space portrait of homoclinic orbits, we include a final sec- 
tion summarizing the generality of zoom-whirl behavior 
and where the separatrix fits in the spectrum of zoom- 
whirl orbits. 



B. Homoclinic Orbits in the Gravitational Wave 
Literature 

For context, we note that homoclinic orbits have ap- 
peared in the gravitational wave literature, although not 
always identified by name. Ref. [l2| analyzed the tran- 
sition for equatorial eccentric Kerr orbits using semi- 
analytic methods. Gravitational wave snapshots and 
semi-analytic estimates of the radiative evolution of or- 
bits near the separatrix appear in [l3|, which also dis- 
cusses the "zoom-whirl" behavior that may be visible 
during an eccentric transition to plunge. The discus- 
sion of separatrices and their role in eccentric transitions 
to plunge is also being discussed for comparable mass 
systems (l5l . H6| . and an eccentric transition to plunge, 
including visible zoom- whirl behavior, has been observed 
in a full numerical relativity simulation of the merger of 
equal mass black holes [I?} ■ 

Homoclinic orbits have also been discussed by name 
in the black hole literature and are not unique to ex- 
treme mass ratio binaries. The distinct imprint on a 
gravitational waveform from the whirl phases or orbits 
near the homoclinic set was discussed in [l8| for both 
Schwarzschild orbits and orbits generated in the Post- 
Newtonian (PN) expansion. A program to identify the 
homoclinic orbits in a higher-order PN expansion is also 
underway (l6| . Ref. [la ] provides a nice summary of 
the interesting phenomenology associated with homo- 
clinic orbits in any dynamical system, and for the case 
of Schwarzschild geodesies formally demonstrates (us- 
ing a somewhat unphysical example) the onset of chaos 3 
around the homoclinic orbits when the system is slightly 
perturbed from the conservative dynamics, a fact that 



Small perturbations to the entire system give rise to structures 
in the phase space, first discussed by Poincare |2(I and usually 
termed "homoclinic tangles", that are quantifiable signatures of 
chaos. 



could be important in the analysis of the transition to 
plunge but which we do not discuss further here. 

While homoclinic orbits are present even in compara- 
ble mass black hole systems described in a PN expansion, 
the complexity of the PN equations of motion makes an- 
alytic results about homoclinic orbits difficult to come by 
for comparable mass systems [? ]. Since most of those 
references compare results against the Kerr equatorial 
case, we restrict our attention here to the fiducial case of 
homoclinic orbits in the equatorial plane of Kerr black 
holes. 



II. ORBITAL PARAMETERS OF HOMOCLINIC 
ORBITS 

The equatorial homoclinic Kerr orbits asymptotically 
approach the same unstable circular orbit in the infinite 
future and past, whirling an infinite number of times as 
they do so. In this section, we provide the afore-promised 
formal definition of a homoclinic orbit and substantiate 
this claim. 



A. Definition of a Homoclinic Orbit 

Formally, a homoclinic orbit approaches the same in- 
variant set in the infinite future as in the infinite past. 
A collection of points S in the phase space of a dynam- 
ical system is an invariant set if orbits that are in the 
set at any time remain in the set for all previous and 
subsequent times. Of course, the set of points in phase 
space traced out by any solution to the equations of mo- 
tion constitutes an invariant set, but useful information 
about global properties of the phase space usually comes 
from identifying invariant sets with some associated re- 
currence property, such as fixed points, periodic orbits, 
or the n-dimensional tori on which bounded quasiperi- 
odic motion in integrable systems unfolds. Henceforth, 
when we refer to an invariant set, we will always mean a 
recurrent invariant set. 

The set of all trajectories that approach S asymptoti- 
cally in the infinite future is a submanifold of the phase 
space, namely the stable manifold of S. Likewise, all tra- 
jectories that approach S asymptotically in the infinite 
past form the unstable manifold of S. A invariant set is 
called hyperbolic if it has both a stable and an unstable 
manifold. 

Now, stable and unstable manifolds of invariant sets 
can sometimes intersect: some individual trajectories 
may approach (possibly different) invariant sets both as 
t — > +00 and as t — > —00. When such a trajectory lies in 
the stable manifold of one invariant set S + and the unsta- 
ble manifold of a different invariant set S-, the trajectory 
is heteroclinic to S + and 5_. If instead the trajectory ap- 
proaches the same invariant set S in the infinite future 
and past, i.e. if it is an intersection of the stable and un- 
stable manifolds of the same set S, then the trajectory 
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is homoclinic to S. 

Identifying the homoclinic orbits in a dynamical sys- 
tem thus amounts to finding the intersections of the sta- 
ble and unstable manifolds of its hyperbolic invariant 
sets. As we will now show, for the system of Kerr equa- 
torial orbits, the only hyperbolic invariant sets with as- 
sociated homoclinic orbits are the energetically bounded, 
unstable circular orbits. Strictly speaking, no relativistic 
orbits are truly recurrent since time itself is a coordinate 
in a relativistic phase space [23|, [28| and all orbits are 
unbounded in their forward motion in time. We will go 
to some trouble in paper II [l4[ to reduce to a 6D phase 
space of spatial coordinates and their conjugate momenta 
in which circular orbits are truly recurrent invariant sets. 



B. Kerr Equations of Motion 

The Kerr metric in Boyer-Lindquist coordinates and 
^eometrized units (G = c = 1) is 



ds 2 = - 1 - — — dt 2 - 



2Mr\ , 9 4Af ar sin 2 9 

?r) dt — 



dtdip 



sin 2 9[r 2 +a 2 



2M a 2 r sin 2 i 



'^d<p 2 , (1) 



A ' 



-dr 2 + Y,d9 2 



of this paper, the first-order geodesic equations are 
Sr = ±VR 



L- 



^ = - (2rE - aL z ) 

A sm a 

(r 2 +a 2 ) 2 E-2arL z . 2 

lit = a hi sm 9 



(3a) 
(3b) 

(3c) 
(3d) 



where an overdot denotes differentiation with respect to 
the particle's (dimensionless) proper time r and 



0(9) 



cos 2 9<a 2 (l- E 2 ) 



L 2 



sm 



R(r) = -(1 - E 2 )r 4 + 2r 3 - [a 2 (l - E 2 ) + L 2 Z ] 
+ 2(aE- L z ) 2 r-QA 



(4) 



(5) 



Modulo initial conditions, we can identify any orbit 
around a black hole of some mass and spin by its con- 
stants of motion E, L z and Q. From this point on, we 
will restrict attention to equatorial orbits. Equatorial 
Kerr orbits have 9 — it/ 2 and 9 = 0. It follows from 
the equations of motion ([3]) that equatorial orbits always 
have Q = and that they remain equatorial and form a 
self-contained set. 

In the following section we will identify the homoclinic 
orbits with the aid of an effective potential picture. 



where M, a denote the central black hole mass and spin 
angular momentum per unit mass, respectively, and 

£ = r 2 + a 2 cos 2 9 

9 • (2) 

A = r — 2Mr + a 2 

Motion along geodesies of ([1]) conserves orbital energy 
E, axial angular momentum L z , the Carter constant Q 
[23 ] , and of course the rest mass ji of the test particle it- 
self. 4 Since there are as many constants of motion as de- 
grees of freedom, the usually second order geodesic equa- 
tions can be integrated to yield a set of 4 first order equa- 
tions of motion for the coordinates [22j ■ Before writing 
them down, we adopt the useful and now common con- 
vention [23| to set both M and fj, equal to 1, an operation 
tantamount to working in units in which the coordinates 
r and t, the proper time r, the spin parameter a and the 
conserved quantities E,L Z ,Q arc all dimensionless. In 
these dimensionless units, which we use in the remainder 



4 E and L z are associated with (-translation and (^-translation 
Killing vectors of the Kerr metric. Q is associated with a Killing 
tensor with a less obvious geometric interpretation. In the weak 
field, Q^Ll + L 2 y . 



C. Effective Potential and the Homoclinic Orbits 

To clarify terms, it is standard parlance to refer to "un- 
stable" circular orbits in the black hole system. Strictly 
speaking, the unstable circular orbits are actually hyper- 
bolic - they possess both a stable and unstable manifold. 
Nonetheless, we continue with this conventional parlance 
to avoid unnecessarily elaborate verbiage and assume the 
reader understands phrasing such as "the stable manifold 
of an unstable circular orbit" . 

The hyperbolic invariant sets in the equatorial Kerr 
system are precisely these unstable circular orbits. Of 
those, the ones that are energetically bound (E < 1) give 
rise to homoclinic orbits. 

Identification of the homoclinic orbits, and indeed in- 
terpretation of the dynamics in general, is easiest with 
an effective potential formulation, which motion around 
spinning black holes admits. However, as we explain 
shortly, the Kerr effective potential has some awkward 
features that make our exposition a bit cumbersome. 
Thus, to ease discussion, we briefly recount the ef- 
fective potential picture for Schwarzschild black holes 
0, HH HE HH and identify the homoclinic orbits for 
that case before extending to Kerr black holes. This sub- 
section amounts to a synopsis of the familiar specifics of 
orbits admitted by the Schwarzschild and Kerr effective 
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potentials, a lengthier accounting of which we include in 
Appendix [A] for the detail conscious reader. 

For Schwarzschild orbits (a = 0), the suggestive form 



1 



-r»-*M=0 
2 2£ 2 



(0) 



of the radial equation (|3a[) becomes the familiar [24j, [2 
1 



-jr 2 + V cH = e cS 



(7) 



describing motion in the one-dimensional effective poten- 
tial 



V eS (r, a = 0) = - 



_ R(r) 



2r 4 



a=0 

2r 2 



2 



(8) 



with effective energy e c s = \E 2 . Note that the asymp- 
totic value of the potential at r = oo is 1/2, so that E = 1 
divides bound from unbound motion. 

An example of such an effective potential for a non- 
spinning black hole with L z = 3.55 is shown in Fig. [2] It 
is simple to read from this figure that the maximum of 
the potential (dV/dr = Q,d 2 V/dr 2 < 0) corresponds to 
an unstable circular orbit and the minimum of the po- 
tential (dV/dr = 0,d 2 V/dr 2 > 0) corresponds to a stable 
circular orbit. Note that the energy of the maximum E u 
is below the asymptotic value E = 1. 

As indicated by the solid line, there is another orbit 
with energy E u but an apastron given by the outer in- 
tersection of the horizontal line of energy E u with V e f{. 
When released from rest at the apastron r a , a test parti- 
cle will roll toward the unstable circular orbit taking an 
infinite amount of time to reach the peak, and likewise 
if time reversed. This orbit is a homoclinic orbit. For 
every bound unstable circular orbit there exists such a 
homoclinic orbit with the same E and L z . 5 Appendix |A1 
shows that these are the only homoclinic orbits. 

For Kerr black holes (a ^ 0), the E and L z depen- 
dences in equation © do not separate as they do in the 
Schwarzschild case. The radial motion can still be cast in 
the form Jjj) as the one-dimensional motion of a particle 
with energy 





£ c ff = 



(9a) 



FIG. 2: R(r) functions and V e «(r) (Eq. @) for 3 
Schwarzschild orbits with L z = 3.55. From bottom to 
top in both diagrams, the corresponding energies are E = 
0.947421, 0.948707 and 0.949993. The first value is for a sta- 
ble circular orbit at r s = 7.679020, the second for an eccentric 
orbit with r p = 6.000593 and r a — 9.656613, and the third 
is for both an unstable circular orbit at r„ = 4.923479 and 
for a homoclinic orbit with r p = r„ and r a = 10.662889. The 
vertical scales have been suppressed for visual clarity. Note 
that R(r) — R'(r) — at the circular orbits. 



moving in a potential 

V cS {r) ee -i?(r)/2£ 2 , 



(9b) 



5 There are no bound orbits with the same (E,L Z ) of unbound, 
unstable circular orbits (i.e. those with E > 1) and therefore 
the unbound circular orbits do not possess homoclinic orbits, as 
elaborated in appendix \K\ 



but unlike when a = 0, the potential depends on both 
E and L z through R(r). V e s is thus a different poten- 
tial for each orbit (i.e. for each (E, L z ) pair) instead 
of a single potential for an entire family of orbits like 
V e s(a = 0). Fig. plots various such functions R(r) in 
the lower panel. As the figure highlights, having the po- 
tential vary under one's feet, so to speak, as the energy 
of the particle changes means that information we could 
previously glean from a single plot of V c ff(a = 0) is now 
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diluted over an infinite number of plots of R(r). Nev- 
ertheless, a bit more effort - expended in appendix [S] - 
shows that even when the black hole spins the unstable 
circular orbits are still the only hyperbolic invariant sets 
and that those with E < 1 give rise to homoclinic orbits. 

Although R(r) changes with E, we can still read qual- 
itative features of the motion effectively from a plot of 
i?(r). To clarify the visual interpretation, Fig.[2]plots the 
R(r) for a = and L z = 3.55 below a plot of the corre- 
sponding Schwarzschild V c g(r). Whereas in the effective 
potential diagram the r values accessible to a particle 
with given energy are those for which V e g is below the 
constant energy line, in the pseudo-potential diagram of 
a given orbit the accessible r values are those for which 
R(r) is above the zero line, reflecting the fact that r in 
Eq. Q is real so that the R(r) under the radicand must 
be non-negative. 

Turning points of the motion, for which V s{r) = s e s = 
0, correspond to single roots of R(r). Circular orbits 
require both V c s = and dV c s/dr = 0, or the equivalent 



R(r) = and R'(r) = 



(10) 



and thus correspond to double roots of R, as Fig. [5] con- 
firms. Simultaneously solving these e quat ions yields ex- 
pressions, originally published in Ref. [271 ], 



E 



L z = ± 



r 3/2 _ 2r l/2 ± , 



r 3/V r 3/2 - 3rV2 ± 2a 

r 2 =p 2ar 1 / 2 + a 2 
r 3/V r 3/2 - 3rV2 ± 2a 



(11a) 
(lib) 



for the energy and angular momentum of circular orbits. 
The top/bottom signs denote progradc/rctrograde. 

Two noteworthy circular orbits deserve mention: the 
innermost stable circular orbit (isco) and the innermost 
bound circular orbit (ibco). As the angular momentum 
decreases, the stable and unstable circular orbits merge 
to a saddle point - the isco. It is the circular orbit for 
which E and \L Z \ arc a minimum 6 271: 



The upshot is that every |L; S co| < \L Z \ < |£ibco| admits 
a bound unstable circular orbit and a corresponding ho- 
moclinic orbit with the same (E, L z ). The apastron of the 
homoclinic orbit with (Ei\, co , ^ibco) is r Q = oo while the 
apastron (and periastron) of the homoclinic orbit with 
(£"ibco) -^ibco) is r a = ^isco- hi other words, the ibco has a 
homoclinic orbit with eccentricity 1 and the isco is a ho- 
moclinic orbit with eccentricity zero. The eccentricities 
of the homoclinic orbits range from 1 down to 0. 



D. Exact expressions for orbital parameters of 
homoclinic orbits 



Above we have described the homoclinic orbits by their 
E or L z . There are other ways to describe the homo- 
clinic orbits. In general, non-circular equatorial Kerr or- 
bits form a two-parameter set, with any particular or- 
bit specified by its energy and angular momentum. For 
bound non-plunging orbits, other pairs of independent 
orbital parameters can also be used, such as the perias- 
tron and apastron (r p ,r a ) or, as is often done, appropri- 
ately defined pseudo-Keplerian parameters (e,p) (eccen- 
tricity and semi-latus rectum, respectively). Homoclinic 
orbits, however, lie in one-to-one correspondence with the 
E u < 1 unstable circular orbits, a one-parameter family 
specified by the radius r u . Homoclinic orbits thus form 
a one-parameter family all of whose orbital parameters 
depend only on the single parameter r u . 

For E and L z this is clearly the case - homoclinic 
orbits have the same energy and angular momentum 
as the circular orbit they asymptotically approach, and 
equations pip determine E and L z once r u is spec- 
ified. Homoclinic orbits also form the separatrix be- 
tween plunging and non-plunging orbits, so they, like 
any bound non-plunging orbit, have well-defined values 
of; 



pi 



" a ,e,p. Simple expressions for those parameters fol- 
low from rewriting the R(r) function, which has a double 
root at r u for homoclinic orbits, as 



R(r) = (E 2 -l)r(r-r u ) 2 (r-r a ) 



(14) 



r isco = 3 + Z 2 T y/(3-Z 1 )(3 + Z 1 +2Z 2 ) ( 12 ) 
Zi = l+ y/l -a 2 [v/TTa"+ VT~a] 
Z 2 = yj'ia 2 + Z\ . 

Since R"{r isco ) = when E = E isco , \L Z \ = \L isco \, the 
isco corresponds to the only possible triple root of R. The 
ibco is the marginally bound E = 1, unstable circular 
orbit H3: 



= 2t« + 2\/T^ 



(13) 



Statements that apply to both prograde and retrograde trajec- 
tories are phrased in terms of \L Z \. 



where r a is the apastron of the homoclinic orbit. Ex- 
panding (fl4|) and equating powers of r with equation (J5]) 
for R(r) yields relations among r u , r a , E and L z . In par- 
ticular, equality of the linear coefficients implies that 



2{aE-L z f 
(1 - E*)rl 



(15) 



Substituting E(r u ) and L z (r u ) from (fTTj) for E and L z 
and simplifying leads to the expression 



2r u (aT 



ru) 2 



4r„ ± 4a 



(16) 



for the apastron of a homoclinic orbit. 

Eq. also furnishes expressions for the e and p of a 
homoclinic orbit in terms of r u . In analogy with Kcplc- 
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rian orbits, the eccentricity 7 and semi-latus rectum of a 
generic orbit are typically denned via 



or equivalcntly 



P 



' a i ' p 

2r a rp 



(17) 

(18) 
(19) 



' a i ' p 

Substituting (fT6|) into (fT8"|) and (fl"9)) with r p = r u yields 
hc _~rl + 6r„ T + 3a 2 



2r u + a 2 



P 



4r u (a =F V^") 2 
ri - 2r u + a 2 



Ref. [l2j derives the implicit relation 

= p 2 (p - 6 - 2e) 2 + a 4 (e - 3) 2 (1 + e) 2 
- 2a 2 p(l + e) [14 + 2e 2 +p(3 - e)] 



(20) 
(21) 

(22) 



that e and p of the homoclinic orbit (referred to there as 
"the separatix") must satisfy, and an equivalent implicit 
expression also appears in [l], [l3| . Acknowledging the 
relationship between the homoclinic orbits and unstable 
circular orbits from the outset furnishes the explicit para- 
metric solutions (|20|) and (|21| to those implicit equations. 



We now know how to specify the equatorial circular 
orbits by a single parameter; either E or L z for instance. 
The unstable circular orbits are a family of hypcrolic sets, 
which means they have stable and unstable manifolds. 
We have also derived the perihelia and apastra of the 
homoclinic orbits as well as the (e,p) as explicit functions 
of r u and spin. 

However, we can do better than this. We can find exact 
solutions for the homoclinic trajectories themselves as a 
function of spin. We will do this now. 



III. EXACT SOLUTIONS FOR EQUATORIAL 
HOMOCLINIC ORBITS 



An exact solution to geodesic motion is a rare commod- 
ity. In this section we very briefly sketch the derivation 
of an exact parameteric solution for homoclinic orbits 
around Kerr black holes of arbitrary spin and refer the 
reader to the acrobatics of appendix f|B]for the detailed 
derivation. 

For any equatorial orbit, the radial motion consists of 
alternating inbound phases [dr/dr < 0) and outbound 
[dr/dr > 0) phases during which r varies monotonically 
with time. Because the equations of motion depend ex- 
plicitly only on r, the radial coordinate parametrizes the 
motion during any single such phase, and other dynami- 
cal variables can be expressed in terms of r. 

Consequently, during an inbound phase, the integrated 
proper time, coordinate time, and azimuth between some 
reference point rg and r are 



r(r) = 
t(r) = 
<p(r) = 



dr— 

dr 

dt/dr 
dr/dr 
dip/dr 
dr/dr 



dr—== 

, r 2 (r 2 + a 2 )E + 2a(aE - L z )r 
dr 

I'D 



dr 



r 2 L z + 2(aE - L z )r 

aTr 



(23) 
(24) 
(25) 



where r and rg are both radial coordinates along the 
same phase (i.e. along a given half-leaf) of the motion. 
Removing the overall minus signs yields the correspond- 
ing expressions for outbound motion. Eqs. (|B1|) - (|B3[) and 
their outbound counterparts are correct for both r < ro 
and r > ro along a single inbound/outbound phase. 



For ordinary eccentric orbits, R(r) has four distinct 
roots and equations (|B1[) - (|F33[) are at best elliptic in- 
tegrals. However, the fact that R factors as in (TH} for 
homoclinic orbits renders the integrals soluble in terms 
of elementary functions. We integrate these equations 
analytically in appendix [B] to give: 
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T 



-(r) 



i(r) 



1 /—, r 2 r a -r 2 , _i T 

--\/r(r a — r H tan W j — tanh w- 

VT^~E^ V V ' (1-^)3/2 V r 7 A r V' 

s y/r(r a - r + 2£^ ^ tan" 1 W- + — tanlT 1 

+ tanh 1 \ - 



-r u r 

r u r a - r 
• a ~r u r 



r a -r 
r a -r + r 



2r_ 

— , tann 
VT^2 



<p(r) = 2 ^ tanh" 1 



r u "Tg-r 



r a - r. 



a _ 
, tann 



r+ r a -r 



Tg-r 



'+ T 



— , tanh 1 i / - 



r a - r 



'a - r_ 



(26a) 
(26b) 

(26c) 



r 



where we have set r = t = ip 
(|T5|) of the homoclinic orbit. 
In the equations above, r 4 



= at r = r a , the apastron 
and r_ represent, respec- 



tively, the outer and inner horizons of the black hole, 
7 = dt/d,T(r u ) and VL U = ^(^«) are the (constant) 
Lorentz factor and azimuthal velocity (£l u > for pro- 
grade orbits, fi u < for retrograde) of the associated 
unstable circular orbit, and E is the energy of the homo- 
clinic orbit (and also of the unstable circular orbit). The 
remaining parameter is 




(27) 



where a prime denotes differentiation with respect to r. 
As derived in paper II in this series [3| , A r is the radial 
stability exponent of the unstable circular orbit. 

Eq ([26) reveals an interesting and unobvious fact about 
the homoclinic orbits. Consider a circular orbit at r = 
r u with energy E and L z and a homoclinic orbit with 
the same energy and angular momentum. Even though 
the homoclinic orbit takes an infinite amount of time to 



asymptote to or away from r u , the total accumulated 
phase difference between the homoclinic orbit and the 
circular orbit over that infinite period is finite. 

To be concrete, consider a prograde homoclinic orbit 
with t = 0, (f = at r = r a , and let the circular orbit at 
r u be at ip = at the same time. Since r varies monotoni- 
cally with t along the homoclinic orbit during its inbound 
phase (as t — ► oo), we can use the r coordinate along 
the homoclinic orbit as a global time parameter via Eq. 
(|26b[) . Since ip along the circular orbit increases linearly 
at a rate lo v = fi u , the phase difference between the cir- 
cular and homoclinic orbits is just the difference between 



cp circ (t(r)) = n u t(r), 



(28) 



and Eq. (|26c|) . By time-reversal symmetry, doubling this 
yields the total phase difference between the circular and 
homoclinic orbits summed over both the inbound and 
outbound phases. Letting t(r) denote time along the 
inbound (f < 0) portion of the homoclinic orbit, the 
resulting phase difference 



A/ c W0) = 2[r (Kr))-V hc (t(r))] 



E 



2 



v / r(r a - r) + 2 



VT^a 2 



3- 2E 2 
1-E 2 

(a — 2Q u r + ) tanh -1 



tan 



1 Ma 



(29) 



r, r„ - r 



r„ - r, r 



(a — 2Q u r_) tanh 1 



r_ r„ - r 



r a — r_ r 



has no divergences, and the limit 

lim A<p hc (t) = lim A</ C (r) (30) 

t — >oc r — >r u 

exists. 

Fig. E] shows how Aip hc depends on r u for various val- 
ues of the black hole spin a. For ease of comparison, 



Aiy9 hc is plotted versus a parameter that varies linearly 
from when r u = r;b C o to 1 when r u = rj sco for a given 
a. The fact that Aip hc ^ mod 2tt for all but a mea- 
sure zero set of homoclinic orbits means that a generic 
equatorial homoclinic orbit asymptotes in the infinite fu- 
ture to a circular orbit out of phase by Atp hc with the 
circular orbit (at the same r u ) to which the homoclinic 
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FIG. 3: The accumulated phase difference between a homocfinic orbit and the unstable circular orbit to which it is doubly 
asymptotic as a function of r„. For a given spin a, the parameter j3 varies linearly from rib co when /3 = to ri sco when (3 = 1. 
Upper: Prograde homoclinic orbits. Lower: Retrograde homoclinic orbits. 



orbit asymptotes in the infinite past. 8 Stated another 
way, if we were to treat all circular orbits at radius r u 
with different phases as distinct, then by adding a con- 
stant and finite phase to any homoclinic orbit, we could 
speak meaningfully about synchronizing it with exactly 
one such circular orbit at t = +oo at r u and with ex- 
actly one circular orbit at t = — oo also at r u but with a 
different phase. 

Except for a measure zero set that accumulate a total 
phase difference Aip hc mod 2ir = relative to a circular 
orbit over their infinite period motion, a homoclinic orbit 
orbit that synchronizes with a given circular orbit at t = 
—oo will be out of phase with that same circular orbit at 
t = +oo by Aip hc . Although a fine detail at this point, 
such phase information could be significant to gravational 
wave templates for the full black hole spectrum. 

We can paint a homoclinic approach to the unstable 
circular orbit using the exact solution of this section. We 
will do so in the context of the special set of periodic 
orbits in the following section. 

IV. HOMOCLINIC LIMIT OF ZOOM- WHIRL 
ORBITS 

Before concluding, we mention another perspective on 
the physical portrait of the homoclinic landmark, and 
that is the connection to zoom- whirl behavior. An associ- 
ation with zoom-whirl behavior has long been suspected, 



This is why we speak about an orbit being homoclinic to some 
invariant set (e.g., the locus of points in phase space with 
r = r u ,p r = 0, ip arbitrary) and not about its being homoclinic 
to a particular orbit (e.g. a particular unstable circular orbit, 
including choice of phase). 



yet also subtley misunderstood. Many practitioners, in- 
cluding the present authors, suspected that zoom-whirl 
behavior was bound to the proximity to the separatrix. 
To the contrary, we found in a previous work [llj . that 
zoom-whirl behavior emerges in the strong-field for any 
eccentricity. Put another way, zoom-whirl behavior is 
demonstrated by orbits that are not in the vicinity of the 
homoclinic orbit as well as by those that are. Still, ho- 
moclinic orbits do have an important significance as the 
infinite whirl limit in the spectrum of zoom- whirl orbits, 
as we now make explicit. 

In Ref. [Il[ , we realized that the spectrum of all black 
hole orbits naturally fall into periodic tables - tables with 
an infinite sequence of entries for a given angular momen- 
tum around a given black hole. Each entry in the periodic 
table is an exactly periodic orbit characterized by a ra- 
tional number that immediately identifies the number of 
zooms, the number of whirls, and the order in which the 
zooms are executed. 

The importance of the periodic tables lies in the ob- 
servation that every bound orbit can be approximated 
to arbitrary precision by a periodic one. Perhaps more 
important to future studies of gravitational waves, every 
bound orbit can be modeled as a slow precession around 
some low-leaf periodic orbit, just as Mercury's orbit can 
be modeled as a precession around an ellipse. 

Although homoclinic orbits are formally aperiodic in 
the sense that they never return to their initial condi- 
tions, they nonetheless mimic periodic orbits, in partic- 
ular they are the infinite whirl limit in the periodic se- 
quence - the final entry in the infinite periodic table pd| . 

To make this connection explicit, consider single-leaf 
orbits like those in Fig. [4j Single-leaf orbits are in one- 
to-one correspondence with the whole numbers; that is, 
the rational associated with each one-leaf orbit counts the 
integer number of whirls. The final entry in this infinite 
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FIG. 4: The progression of the 1-leaf periodic orbits through 1,2,3, 4... oo whirls. The orbits shown are prograde orbits for 
a = 0.5 and L z — 3.158540, or the average of Li sco and iibco- Note that the whirls beyond the second whirl are too closely 
packed in r to distinguish visually in the plot. 



list is the orbit that executes an infinite number of whirls 
and therefore never actually reaches the end of its first 
radial cycle. That orbit is of course the homoclinic orbit. 



More specifically, if we write down a sequence of en- 
ergies E w , eccentrities e w , apastra r aw , periastra r p , 
actions J rw , G ^ c - f° r the constant L z set of one-leaf or- 
bits, then all of these sequences converge in the w —> oo 
limit to the values for the homoclinic orbit with the same 
L z (for a black hole of a given spin). 



As already noted, for a given a and L Zl the homoclinic 
orbits form the separatrix in the phase space between 
orbits that are energetically bound and those that are 
not. As discussed in detail in [ll|, [23|, [28[, orbits that 
are bound in the phase space turn out to lie on surfaces 
homeomorphic to 2-dimensional tori (3-tori for generic 
nonequatorial orbits), and only these bound orbits have 
an associated set of fundamental frequencies in terms of 
which orbit functionals can be Fourier expanded [|| . The 
homoclinic orbit of a given L is also therefore the sepa- 
ratrix between the regions of phase space inhabited by 
these quasiperiodic orbits and those that are fully aperi- 
odic. Since all quasiperiodic orbits can be approximated 
by the periodic set, the homoclinic orbit is the divide be- 
tween the domain of influence of the periodic set with its 
correspondence to the rationals and aperiodic orbits that 
merge or escape. 



V. CONCLUSIONS 

Homoclinic orbits offer the kind of crucial signpost that 
demarcates physically distinct regions of the conservative 
and inspiral dynamics: bound from plunging, whirling 
from not-whirling, smooth from chaotic. They thereby 
define salient details of black hole dynamics and we have 
spent time in this article deriving an exact parameteric 
solution for homoclinic motion that we hope will prove 
of use to others in the field. 

Physically, we have shown that homoclinic trajectories 
are an infinite whirl limit of the zoom-whirl orbits. Even 
remembering that zoom- whirl behavior is generic and not 
exotic in the strong- field [ll[ , the homoclinic orbits them- 
selves are a special and sparse subset. Nonetheless, every 
inspiraling orbit must transit through a homoclinic orbit 
on the transition to plunge. The isco, which is the exit to 
plunge for quasi-circular inspiral, is itself a homoclinic or- 
bit with eccentricity zero. The homoclinic family ranges 
in eccentricity from zero (the isco) all the way up to 1 
(homoclinic to the ibco). All orbits, except those ex- 
ceptionally well-approximated as quasi-circular, will roll 
through another member of the homoclinic family on the 
transition to plunge. 
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APPENDIX A: THE EFFECTIVE POTENTIAL 

1. a = 

Although the key features of the Schwarzschild geome- 
try are recognizable at a glance to anyone familiar with an 
effective potential formulation, the Kerr case is less visu- 
ally informative. In order to ground the details of the cir- 
cular and homoclinic orbits, we include the Schwarzschild 
treatment in detail. 



^effM 




FIG. 5: The Schwarzschild effective potential drawn in solid 
lines as a function of radial coordinate r for L z = 3, 3.8 and 
4.4, from bottom to top. The upper and lower dashed lines 
represent the borderline potentials for L z = Libco and L 2 = 
iisco, respectively. 

The value of L z fixes the form of the potential, as 
shown in Fig. El Two critical values of L z define regimes 
in which the potential exhibits different qualitative fea- 
tures: the angular momentum Lj SCO of the innermost sta- 
ble circular orbit (isco) , associated with the saddle point 
in the lower dashed potential, and the angular momen- 
tum Libco of the innermost bound circular orbit (ibco), 
the circular orbit with E = 1. 

The effective potential picture allows us to determine 
the hyperbolic invariant sets at a glance: 

• No invariant set exists with L z < Li sco , since ev- 
ery such orbit plunges and thus fails the recurrence 
test. 



• When L z > Li sco , the potential admits one stable 
circular orbit with r = r s , E = E s (minimum of 
V c g(a = 0)) and one unstable circular orbit with 
r = r u , E = E u (maximum of V e g(a = 0)). 

• The value Libco further distinguishes the two sub- 
cases seen in Fig. [6j from which we see that orbits 
with E E u never approach an invariant set. In- 
stead, as t — » ±oo, every such orbit (a) oscillates 
between two turning points in the potential well, 
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L z > L ibc 




L isco < L z < L ihco 




In contrast, orbits with i2 = _E U do asymptotically 
approach an invariant set, namely the unstable circular 
orbit. Consider the upper panel in Fig.[6j for which L z > 
iibco- There is an E = E u orbit that approaches r u 
as t — > — oo and plunges as r, — > +00 and another that 
plunges as t — > —00 and approaches r u as i — > +00, both 
represented by line (f) in the figure. These two orbits 
are distinct, just as the E = E u orbit that escapes as 
t — > — 00 and approaches r u as t — > +00 is distinct from 
its time-reversed counterpart (both represented by (g)). 
So while they define stable and unstable manifolds for 
the circular orbit shown, these E = E u orbits are non- 
intersecting (share no initial conditions (r, f)) and thus 
are not homoclinic to the circular orbit. 

However, when Li sco < L z < L[\ JCO , as in the lower 
panel of Fig. the E — E u orbit (h) has a turning point 
and thus approaches r u at both t — > ±00. Parts of the 
stable and unstable manifolds of this unstable circular 
orbit intersect (in fact, they completely coincide), and 
these orbits are therefore homoclinic to the circular orbit. 

We thus conclude that since they are the only recur- 
rent orbits that are approached by any other orbit in the 
infinite future or past, the unstable circular orbits are 
the only hyperbolic invariant sets. Furthermore, those 
unstable circular orbits with E < 1 (i; SC o < L z < Ln, co ) 
have associated homoclinic orbits with the same angu- 
lar momentum and energy, or more specifically, a family 
of such orbits differing from one another by an overall 
translation in (p. References [l!| and [l8[ make similar 
arguments for the Schwarzschild case. 



2. a 



FIG. 6: Representative plots of the Schwarschild V e g(a — 
0) for a value of L z > Libco (top) and Li sco < L z < Lvaco 
(bottom). As explained in the text, the horizontal lines define 
the energies (for the fixed L z of each V e s(a = 0)) of orbits 
that (a) oscillate, (b) plunge, (c) escape to r = 00, or (d) 
both plunge and access r = 00, as t — » ±00. The lines (f) and 

(g) tangent to Vesip, — 0) at r = r u represent E = E u orbits 
that asymptotically approach r u at t = +00 or —00 (the ones 
marked (f) also plunge). In the lower figure, the E — E u orbit 

(h) on the right also has a turning point, so it approaches r u 
at both t — ±00 and is homoclinic to the unstable circular 
orbit. 

(b) plunges, 9 (c) escapes, or (d) escapes as t — > —00 
and plunges as t — > +00 (or vice versa) . 

None of the above asymptote to an invariant set. 



Eq. d 3d I ) implies that orbits plunge (reach the horizon) after a 
finite amount of proper time r but an infinite amount of coordi- 
nate time t. 



Our argument will focus on the roots of the quartic R, 
which we can rewrite as 

i?,(r) = (E 2 - l)r(r - n){r - r 2 ){r - r 3 ) . (Al) 

For ease of notation, we adopt the conventions that, from 
left to right in (|A1|) . real roots appear before complex 
roots and the nonzero real roots appear in ascending or- 
der n < r-2 < r^. Additionally, 

R'(r = 0) > (for aE ^ L z ) , (A2) 

so R is negative just to the left and positive just to the 
right of the root at r — 0. Since complex roots occur in 
conjugate pairs, the zero root implies that at least one of 
the three remaining roots is real. 

The non-negativity of f 2 implies that motion is only 
possible where R(r) > 0. Turning points of the mo- 
tion, for which V c s{r) = £ c ff = 0, correspond to single 
roots of R(r). Circular orbits require both V e e = and 
dV e s/dr = 0, or the equivalent 

R(r) = and R'(r) = , (A3) 
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and thus correspond to double roots of R, as Fig. [^con- 
firms. Simultaneously solving these equations yields ex- 
pressions [13 



r 3/2 _ 2r l/2 ± , 



r 3 /Vr 3 / 2 - 3rV2 ± 2a 



± 



" 2 T 2ar 1 /2 



r 3/4^/ r 3/2 _ 3r i/ 2 ± 2a 



(A4a) 
(A4b) 



for the energy and angular momentum of circular 
orbits, where the top/bottom signs apply to pro- 
grade/retrograde orbits. These functions, plotted for a 
sample of a values in Fig. W\ simultaneous minima (max- 
ima for retrograde L z ) at [27| 



r isco = 3 + Z 2 T y/(3-Z 1 )(3 + Z 1 +2Z 2 ) (A5) 
Zx = l+ Vl-a 2 [\/T+~^ + \/T~^\ 
Z 2 = ^3a 2 + Z 2 . 

Since R"(r isco ) = when £ = E isco , \L Z \ = \L isco \, the 
isco corresponds to the only possible triple root of R. 
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FIG. 7: E and L z as functions of the radius r of prograde (left 
panels) and retrograde (right panels) circular orbits. The spin 
parameter for the curves are a = (solid curve), and then in 
order of increasing distance from the solid curves, a = 0.8, 0.9 
and 0.995. For a given a, E and L z have simultaneous minima 

at T — 7"i sco . 



Furthermore, since 

d 2 V eS 
dr 2 



R =o, 
R'=o 



R"{r) 
2T? 



(A6) 



(A7) 



R" also determines the stability of circular orbits, with 

r < ri SCO ==> R"{ r ) > ==> unstable 
r > r iBCO R"(r) < stable 

for a given \L Z \ > |Li sco |- Also, paralleling the a = 
case, E(r) in (|A4a[) increases monotonically for r > r- lsco 
and approaches 1 as r — > oo, so that stable circular orbits 
always have Ei sco < E < 1. Unstable circular orbits, on 
the other hand, can have any E > Ei SCO , and the circular 
orbit with E = 1 occurs at 1271 



2=Fa + 2VlT^ 



(A8) 



We now show that every non-circular Kerr equato- 
rial orbit falls into one of the same categories listed for 
Schwarzschild orbits in Fig. [S] Recall that each R(r) 
plot represents only those orbits with the same E and L z 
and that motion is only possible in regions where R > 0. 
When E > 1, R(r) -> +oo at both r -> ±oo, and (|A"2]l 
implies that i? has a negative root. There arc thus three 
possibilites for the number and type of positive roots: 

• No positive roots, in which case all positive r are 
accessible, and R represents a single type (d) orbit. 

• Two positive roots r 2 < r$, resulting in a type (b) 
orbit (0 < r < r 2 ) and a type (c) orbit (r > r$ ). 

• One positive double root r 2 = r% = r u with 
R"(r u ) > 0, in which case R represents an E u > 1 
circular orbit at r = r u plus orbits of type (f) 
(0 < r < r u ) and (g) (r > r u ) that asymptot- 
ically approach r u in either the infinite future or 
past (but not both). 

The E = 1 case is the same as above but without the 
negative root (since R is only cubic when E = 1). We 
thus conclude as in the a = case that the invariant sets 
with E > 1 are unstable circular orbits but that they do 
not have orbits homoclinic to them. 

When E < 1, R — > oo at r - ► ±oo. Eq. (|A2[) requires 
that there be at least one positive root and only an even 
number of negative roots, leaving four possibilites for the 
number and type of positive roots: 

• Just the one root r\, resulting in a type (b) orbit 
(0 < r < n) 

• Three total positive roots r\ < r 2 < r^, resulting 
in a type (b) orbit (0 < r < r\) and an oscillatory 



bound orbit of type (a) (r 2 



1 p 



< T < ^3 



1 a) 



• One single positive root T\ and one double root 
r 2 = r s > r\ with R"(r s ) < 0, denoting a type (b) 
orbit (0 < r < r\) and a stable circular orbit of 
radius r s 

• One single root r 2 and one double root at T\ = 
r u < r 2 with R"(r u ) > 0, so that R{r) features 
an unstable circular orbit with E u < 1 at r„, a 
type (b) orbit (0 < r < n), and a type (h) orbit 
(r u < r < r% = r a ) that approaches r u as t — > ±oo, 
i.e. an orbit homoclinic to r u . 
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As before, we conclude that the invariant sets with homo- APPENDIX B: DERIVATION OF EQUATORIAL 
clinic orbits are the unstable circular orbits with E u < 1, HOMOCLINIC ORBITS 

one of which exists for every ribco < r u < ri sco . 

We now find an exact solution by analytically integrat- 
ing the following equations: 

^-jOs-jM (B1) 



f(r) = - / dr dt/dT 



dr/dr 

, r 2 (r 2 + a 2 )E + 2a(aE - L z ) 
dr == 



ip(r) = — j dr 



dr/dr 

, r 2 L + 2(aE- L z )r 

dr == — 

Av^R 



(B2) 



(B3) 



where r and ro are both radial coordinates along the 
same phase (i.e. along a given half-leaf) of the motion. 
Removing the overall minus signs yields the correspond- 
ing expressions for outbound motion. 



1. Integral equations for homoclinic orbits 



Some algebraic manipulation of the denominators of 
the integrals (|B1[) - (|B3[) renders them more suitable for 
evaluation. We begin with R(r), which for equatorial 
orbits has its smallest root at r = (since Q = 0). For 
homoclinic orbits specifically, the remaining roots of R 
are a double root at r u (= r p ) and a simple root at r a . 
R(r) therefore factors into 

R(r) = (E 2 -l)(r-r u ) 2 r{r-r a ) (B4) 
= (l-E 2 ){r-r u ) 2 r{r a ~r), (B5) 

where we've written R in the second form so that the 
product of the r-dependent terms is manifestly positive 
for r p < r < r a . The square root in the denominators 
thus becomes 



y/W) = yf\ - E 2 (r - r u )^Jr{r a ~ r) , (B6) 



where we have replaced \/ {r — r u ) 2 with (r — r u ) since 
r > r u over the entire orbit. A also factors into 

A= (r-r+)(r-r_) , (B7) 



where r + = 1 + \J\ — a 2 and r _ = 1 — y/l — a 2 are the 
outer and inner horizons, respectively, of the central black 
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hole. The integrals (|B1[) - (|B3|) are therefore 
1 



r(r) 



t(r) 



dr 



dr 



ip(r) = - 
di 



Vi 


-E 2 




C 




1 


Vi 


-£ 2 


) 




(r- 


-r + )( 




1 



(r - r u )yjr{r a - r) 



r 2 (r 2 + a 2 )E + 2a(aE - L z )r 



r 2 L z + 2(aE - L z )r 



(B8) 



(B9) 



(BIO) 



(r - r+)(r - r_)(r - r u )y/r{r a - r) 



2. Change of variable 

We can express the integrals above more compactly as 



r(r) 
t(r) 



1 



Vi- 


£ 2 


1 




Vi- 


£ 2 






1 





/l 



p(r) = [L 2 J 3 + 2(aS - L Z )I 4 



vnvE 2 



(BH) 
(B12) 
(B13) 



where 



/a = ~ dr 



dr 



(r 


- r u )y/r(r a - r) 

( 




r 4 






(r - r+)(r 


- r_)(r - r u )^/r(r a - r) 

( 








r 2 


(r - r+)(r 


- r_)(r - r u )^/r(r a - r) 

1 




r 



(B14b) 
(B14c) 



(r - r+)(r - r_)(r - r u )yjr{r a - r) 



(B14d) 



We evaluate each of the integrals in (|B14j) in closed form 
by the same procedure. First, we bring one (positive 
definite) r from each numerator under a radical as an r 2 : 



h 
h 
h 

h 



dr 
dr 



r 



(r - r u ) 

^.3 



dr 



dr 



r a - r (r - r + )(r - r_)(r - r u ) 
r 



r 



r a - r (r - r + )(r - r_){r - r u ) 



1 



r a — r(r- r+)(r - r_)(r - r u ) 
Next, the change of variable 



(B15a) 
(B15b) 

(B15c) 
(B15d) 



ir a — r r 
u = \l , r 



—dr 



= du ■ 



u 2 + 1 
2r a 



(B16) 



r a — r (u 2 + l) 2 

turns the integrals (|B15j) into 



h = 2r 2 / du 



Io = 2r 4 / du 



1 



y 2 {r a - r u y) 



(B17a) 



«„ y 2 {r a -r + y){r a -r_y){r a -r u y) 



(B17b) 



7 3 = 2rf / cfc 



1 



(r a - r + y) (r a - r_y) (r a - r„y) 



(B17c) 



(B14a) J, =2r„ / du 



(r a - r + y) (r a - r_y) (r a - r u y) 



(B17d) 



where we've written y = u 2 + 1 as a shorthand. 



3. Partial fraction decomposition 



Each integrand in (|B17[) is now a product of factors 
linear in y and splits up via a standard partial fraction 
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decomposition. We get 



where 



h = 2 du 



Aii 
y 2 



A12 

V 



T a - r u y 



(B18a) 



h = 2 / du 
' itt) 



A21 _j_ A 2 2 
,,2 



.4 



2:i 



V r a - r u y 



A 
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r a -r+y r a -r_y 




(B18b) 



(B18c) 



(B18d) 



An=r a =A 2 i A12 — r u A 2 2 = r u +r + + r- A X3 = r u 

A23 = rlA± 3 A33 = r u Ai 3 



A 2i = r A + A AA A 34 = r + A M 
A 35 = r_A 



+ 

A 25 = r 3 A 45 



45 



A43 = 
A iA = 
A 45 = 







- r+) (r u 


-r_) 




~r+ 




(r + 


- ?'_) (r u 


-r + ) 




r_ 




(r + 


- r_ ) (r u 


-r_) 



(B19) 



We are left with five different integrals to calculate. Re- 
calling that y = u 2 + 1, those antidcrivatives evaluate to 



Ii = 



du = f du- 



(u 2 + iy 



2 \l + u 2 



[If 1 

I2 = 1 du — = du —= 

J y J u 2 + i 

= tan -1 u 



T-i — I du 



1 



r a - r u y 
1 



du ■ 



(B20a) 



(B20b) 



r a - r u - r u u" 



fuy/^a Til 

du 

r a - r + y 

1 



tanh 1 < u 



a 1 u 



du ■ 



(r a - r + ) - r + u 2 



; t anh 



X5 = du 



1 



r a - r_y 
1 



du ■ 



{fa — r _) — r _u 2 



tanh 1 i u 



(B20c) 



(B20d) 



(B20e) 
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The integrals (|B18[) are therefore 



where the functions fj(r) are 



h 


= 2 (A n Ii - 


h A12J2 ■+ 


■A13I3) 


(B21a) 


h 


= 2 (A 21 J 1 - 


f- A22X2 4 
f A25I5) 


- A23T3 


(B21b) 


h 


= 2 (A 33 l 3 - 


1- A34I4 4 


-A 35 I 5 ) 


(B21c) 


h 


= 2 (A 43 T 3 - 


h A44X4 4 


- A45T5 ) 


(B21d) 



Combining equations ([BTT|l - |(BT3|) . (|B19]) . |B20|) and 



(|B21[) , and recalling that u = ^/ (r a — r)/r, the dynami- 
cal variables r, i and </j become 



r(r) 
t(r) 



1 


Vi- 


£ 2 


1 




Vi- 


E 2 


1 





y/l-E 2 



3=1 

3=1 



(B22) 
(B23) 
(B24) 



Mr- 
Mr 

/s(r 

/t(r 

/s(r 



= \Jr (r a - r) 



tan 



tanh 



tanh 



tanh 





r v 








r 


r a 




r„ 




r 






''"+ 








r 


r a 




''■+ 




r 






?•_ 








r 


r„ 




r_ 




r 





(B25a) 
(B25b) 

(B25c) 

(B25d) 

(B25e) 



3=3 



and the corresponding coefficients are 



c{ r) = 1 Cf' = (r„ + 2r u ) 

cf = E 



d T) = 2. 



E (r a + 2(r u + 2)) 



C, 



a) 



(/) 



a 'it 



G 



(¥>) 



X 



a 'it 



r£ (rj + a 2 ) £? + 2a (a£ - L.) r u r 2 u L z + 2(aE~ L z ) r u 



(r u - r + ) (r„ - r_)r 2 



-4r 4 



(r„ - r + ) (r n - r_)r 2 
-2r + 



(B26) 



C 



(/) 



2£V + - aL z 



yjr + (Fa - r+) (r u - r+) 
4r_ 



2a£; - £ 7 r_ 



V r + ( r a - » + ) - r+) 
2r_ 



2Er_ - aL 2 



V r - ( r a - r_) (r„ - r_) 



2a£; - L z r + 



V r - ( r a - r_) (r u - r_) 



In the expressions above, we have used the facts that 
(since r ± are roots of A) r 2 +a 2 = 2r ± and that r + +r_ = 
2. 

Note that since all of the fj(r) vanish at r = r a , our 
expressions (|B22|) - (|B24|> implicitly assume the natural 
choice of time and azimuthal origins, namely at apastron. 
After all, since they are single-leaf orbits with formally 
infinite radial periods, homoclinic orbits have only 1 out- 
bound and 1 inbound phase each and transit through 
apastron only once. We now make that choice explicit. 



From here on, all expressions assume that 

r(r a ) = t(r„) = <p{r a ) = , (B27) 

along homoclinic orbits, so that r and t are posi- 
tive/negative along the inbound/outbound branch, while 
ip is positive/negative along the inbound/outbound 
branch for prograde orbits (increasing tp) and nega- 
tive/positive along the inbound/outbound branch for ret- 
rograde orbits (decreasing ip). 
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4. Simplification of coefficients 



The coefficient ci^ is therefore 



The task now is to render the coefficients in a more 
meaningful form. The Ci's are already simple. To sim- 
plify the C2's, we expand the factored form (|B5[) of 
R(r)to 

R(r) = {1-E 2 )x 

{-r 4 + (2r u + r a )r 3 - r u (2r a + r u )r 2 + r 2 u r a r) 

(B28) 

Comparing to (|5|) (with Q = 0) and equating coefficients 
of corresponding powers of r, we see that for equatorial 
homoclinic orbits, 



C. 



(r) 



E 2 



(B35) 



C { 3 t] and are each ci^ times another factor. Com- 



paring to (|B1[) and (|B3|) . however, we can identify these 
extra factors as the dt/dr and dip/dr, respectively, of the 
unstable circular orbit associated with the homoclinic or- 
bit. That allows us to write 



Thus, the C2's are 



r {r) _ * 



1-E 2 



a 



(*) _ 



E 



= 2E 



1-E 2 
3-2E 2 
1-E 2 



(B29) 

(B30) 
(B31) 



For the CVs, notice that 

R"{r) = (1 - E 2 ) x 

{-12r 2 + 6(2r u + r a )r - 2r u {2r a + r u )} 

(B32) 

=^ R"(r u ) = (l-E 2 u )2r u (r a -r u ) . (B33) 

Inserting this into the expression for the proper time sta- 
bility exponent jX r of the unstable circular orbit associ- 
ated with the homoclinic orbit yields 



'R"(r u ) 
2E?, 



'2r„(r a 



2rt 



(B34) 



a 



(0 



■dt 



^-^E 2 ^- (r„) = f v/T^ , (B36) 

where A r , as we will show in paper II |14| . refers to the 
stability exponent governing the evolution with respect 
to coordinate time t of small perturbations to the circu- 
lar orbit (we've used here the fact that X r dt = jXrdr). 
Likewise, 



a 



(¥>) 



— Jl^E^lr ) - - — (r Ar ) 

7 A r V dT [u> X r dt (u 'dT [u) 



Ar Clt 



(B37) 



where Q u = ^{r u ). 

Simplifying the CVs take a little more work. As men- 
tioned in § HH the energy and angular momentum of 
the homoclinic orbit are the same as those of the un- 
stable circular orbit at r u . Recalling the expressions 
(|A4bp for circular orbits (top/bottom signs are for pro- 
grade/retrograde orbits) from [27| . we can rewrite the 
numerator of the second factor in cjjp as 
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2Er + — aL z 



2r + (rl /2 - 2r]j 2 ± a) T a (V 2 qp 2ari /2 + a 2 

3/4 / 3/2 ~ 1/2 i o 



3/2 . 1/2 . 

2r + r u — 4r + 7V 4 


-2^r.y 2 




r 2 - 2r + + a 2 












/ 3/2 „ 1/2 , „ 

/ tu - 3r u ± 2a 



o 1/2 

2r + r u ' 



"„ - 2 + r_ =f a (r u - r + ) (r„ + r + ) 



(B38) 



( r u -r+) 



3/4 4 / 3/2 o 1/2 i n 

r u V rV - 3r u ' ± 2a 



1 /2 

2r + 7V To(r„ + r + ; 



r 



where we've used the facts that r± are roots of A and that 
r + r_ = a 2 . Analogously, the numerator of the second 
factor in Cg becomes 



2Er_ - aL z 



(r„ - r_ ) 2r_ 2 =F a (r u + r_ ) 



3/4 / 3/2 „ 1/2 , „ 

r„ V r u ' - 3r u ' ± 2a 



leaving the coefficients C± \ as 

-4r . 



d* 5 



/2 -3ri /2 ±2a 



3/2 / 3 
n 1/2 / \ 

2r + r„ To(r„ + r + ) 



\/r~+ ( r a - r + ) 



a 



(t) 



4r_ 



x 



3/2 / 
1/2 



3/2 „ 1/2 , o 

r« - 3r„ ± 2a 



(B39) 



(B40) 



(B41) 



2r_r„ =pa(r u +r_) 
x — - 

y/r- (r a - r_) 

For what follows, it will be useful to look at the signs 
of the numerators 



2r + r 1 J 2 Ta{r u +r + ) , for C* 4 
2r_r 1 J 2 Ta(r u +r_) , for C* 5 



(B42) 
(B43) 



In the retrograde case (bottom sign), each is the sum of 
two non-negative terms and thus strictly non-negative. 
In the prograde case (top sign), we can whether the sign 
depends on the values of r u and a by treating each of 
(|B42|) and (|B43[) as quadratic function of the variable 

y = rl/ 2 . Specifically, those functions will be negative 
when 



ay 2 — 2r + y + ar + > , for C4 
ay 2 - 2r_y + ar_ > , for C 5 



(B44) 
(B45) 



Since the expressions above have positive quadratic co- 
efficients, the inequalities (|B44[) and (|B45[) are satisfied 
when 



y < ^ (l - v/l-r_) or y > r -± (l + 



for Ca and 



(B46) 



V < 



or y > 



1 + 



(B47) 

for C5. In the case of (|B47|) . the radicand 1 — r + = 
— yl — a 2 is strictly negative 10 and the roots are com- 



10 For a = 1, the radicand is 0, not negative. However, in this 
scenario, r— = 1 and the quadratic expression in ||B45|I has a 
double root at y = 1 r u = 1. Since r u > 1 for all a, 

then even in the a = 1 case, the quadratic in I IB45D will be non- 
negative. Of course, the a = 1 case for any analysis of orbital 
motion must be handled carefully since the r coordinate values 
of the inner and outer horizons, the itco, the ibco and the isco 
are all unphysically degenerate in the maximal spin case. 
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plex. The quadratic expression is therefore always posi- 
tive, and (|B43|) is always negative. 
For (|B46j) . we note that since 



A simple plot (not included here) shows that 



< 1 for < a < 1 



(B48) 



the lower root is subhorizon and thus irrelevant (because 
T u > t+ ) ■ So what we must check is whether we can ever 
have 



< 



for < a < 1 



(B51) 



2 r + 

y > -j 



(B49) 



In fact, (|B49|) is never satisfied for prograde orbits. To 
see why, recall [23] that for prograde equatorial orbits, 



3 + Z 2 -[(3- Z 1 )(3 + Z 1 +2Z 2 )} 



1/2 



Zx = 1 + (1 - a 2 ) 1/3 [(1 + a) i/a + (1 - a) 

.2 , ^2\l/2 



1/3 



Z 2 = (3a 2 + Z 2 ] 



(B50) 



Since r u < rj SCO for all (eccentric) homoclinic orbits, 
(|B49j) is never satisfied, and (|B42|) is always positive. 



The upshot is that we can write the C^l 's so that every 
factor outside a radical is positive. Specifically, 



(t) 



G 



(/) 



r> 1/2 / \ 

2r + 7V =F a (r„ + r + J 



\jr + [2-{2r u +r + ){l-E 2 )] r 3 / 2 ( 



3/2 „ 1/2 , o 

r« - 3r„ ± La 



Vl-E 2 



— 2r_rv ± a (r M + r_) 



[2 - (2r u + r_) (1 - E 2 )] r 3 / 2 (r 3/2 - 3ri /2 ± 2a) 



(B52) 



(B53) 



where we have used equation (|B29j) to rewrite the factors Since the numerators of the factors on the second lines are 



(fa ~ r ± ) m the radicands of the denominators as 
_ 2-(2r M + r ± )(l-£ 2 ) 



(B54) 



now manifestly positive, they can be brought under the 
radical sign without having to worry about stray factors 
of —1. We are left with 



Gf 



4r i 



G, 



(/) 



4r_ 



Vl-E 2 





"r, 1/2 / n" 

2r + r u T a (r u + r + ) 


2 


r + [2 - (2r„ + r+) (1 - E 2 )] r 3 J 2 (r 3/2 - 3r!/ 2 ± 2a) 




-2r_rV 2 ± a(r u + r_) 


- 2 



\ r_ [2 - (2r„ + r_) (1 - £ 2 )] r 3/2 (r 3 



3/2 -3ry 2 ±2a 



(B55) 



(B56) 



Finally, each of the large radicands in (|B55[) . (|B56|) is 1. in each denominator as 
To see this, we use equation (|A4a[) to rewrite the 1 — E 2 

l-E-= r "~ 4r ; ±W ~ a \ (B57) 
rT (rT - Zrl /2 ± 2a) 
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so that distributing the r 3 / 2 (r 3 / 2 — Sri/ 2 ± 2a ) in the 



denominators of the radicands leaves them in the form 
Wj 2 (rT - 3rV 2 ± 2a 



(2r u + r + ) (r 2 u - 4r„ ± Aar]J 2 - a 2 ) 



for C, 



(t) 



(B58) 



2r 3 J 2 [r 3 J 2 - 3rl/ 2 ± 2a 
-(2r„ + r_)(r£ 



4r„ ± ±ar l J 2 



for a 



(t) 



(B59) 



Multiplying out the numerators and denominators and 
grouping them by powers of r u then shows that they are 
identical, for both prograde and retrograde orbits. 

The final expressions for the coefficients C^l are com- 
pact. Noting that r + — r_ = 2yl — a 2 , those expressions 



a 



2 I±= x v/T~^ 



t) 



VT^ a 

2r_ 



\/l-E 2 



(B60) 
(B61) 



To get the corresponding ip coefficients, note that 
1 



2aE — L z r^ 



2aEr ± — L z r^r ± 



(B62) 



Thus, 



r<(v) _ 

° 4 ~2~ La ~ 
r M _ 1 a r (t) _ 



1 a 



it) 



\fl-E 2 
\/l-E 2 



(B63) 
(B64) 



To summarize, once simplified, the coefficients in (|B26[) 
become 



Cf > = 1 



cf 5 



1-E 2 



a 



(r) 



2 



cf = E 



C?> = 2E 



C 



(i) 



2r+ 



a 



'A, 



y/i 


-E 2 


'3- 


2E 2 \ 




-E 2 J 


y/l 


-E 2 



VT^a- 



:^\-E 2 a 



c 



2r_ 



s/T^a 2 



Vi-e 2 a 



(¥>) 

5 



VT^ 2 



(B65) 
V^-E 2 

\Jl-E 2 



5. Analytic expressions for homoclinic orbits 

We can now put everything together from the prior 
subsections. Looking back at equations (|B22p - (|B24|) 
and substituting from (|B25P and (|B65[) . we arrive at the 
final expressions for all the dynamical variables: 



— [2Er ± - aL z ] 
~r± 



r(r) 
f(r) 



1 



VT^e 2 

E 

2r+ 



\J r(r a - r) + 



(1-E 2 ) 



2\3/2 



tan 



1 Ha 



y/r{r a -r) + 2E-^- ^4 



(1-E 2 ) 



3/2 



r 7A r 

2 



tanh 



r u r„ - r 



r„, - r u r 



Vl - a 2 



tanh 1 



r + r a — r 2r_ 



Lp(r) = 2—— tanh W 

Ar \ T a — T u T 



r a ~r + r \/l - a 2 
a 



r X r 
tanh -1 



tanh 



r„ r„ - r 



r„, - r,. r 



r_ r„ - r 



r a -r_ r 



VT^a 2 



: tanh 1 



r , r„ - r 



r a — r+ r 



tanh 1 



r_ r„ - r 



r a -r_ r 



(B66) 



(B67) 



(B68) 
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